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Abstract Clusters of galaxies are self-gravitating systems of mass ~ 
and size ~ 1 — 3h~ l Mpc. Their mass budget consists of dark matter (~ 80%, on average), 
hot diffuse intracluster plasma (< 20%) and a small fraction of stars, dust, and cold gas, 
mostly locked in galaxies. In most clusters, scaling relations between their properties, like 
mass, galaxy velocity dispersion, X-ray luminosity and temperature, testify that the cluster 
components are in approximate dynamical equilibrium within the cluster gravitational po- 
tential well. However, spatially inhomogeneous thermal and non-thermal emission of the 
intracluster medium (ICM), observed in some clusters in the X-ray and radio bands, and the 
kinematic and morphological segregation of galaxies are a signature of non-gravitational 
processes, ongoing cluster merging and interactions. Both the fraction of clusters with these 
features, and the correlation between the dynamical and morphological properties of irreg- 
ular clusters and the surrounding large-scale structure increase with redshift. 

In the current bottom-up scenario for the formation of cosmic structure, where tiny fluc- 
tuations of the otherwise homogeneous primordial density field are amplified by gravity, 
clusters are the most massive nodes of the filamentary large-scale structure of the cosmic 
web and form by anisotropic and episodic accretion of mass, in agreement with most of the 
observational evidence. In this model of the universe dominated by cold dark matter, at the 
present time most baryons are expected to be in a diffuse component rather than in stars and 
galaxies; moreover, ~ 50% of this diffuse component has temperature ~ 0.01 — 1 keV and 
permeates the filamentary distribution of the dark matter. The temperature of this Warm-Hot 
Intergalactic Medium (WHIM) increases with the local density and its search in the outer re- 
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gions of clusters and lower density regions has been the quest of much recent observational 
effort. 

Over the last thirty years, an impressive coherent picture of the formation and evolution 
of cosmic structures has emerged from the intense interplay between observations, theory 
and numerical experiments. Future efforts will continue to test whether this picture keeps 
being valid, needs corrections or suffers dramatic failures in its predictive power. 

Keywords Galaxies: clusters: general ■ Intergalactic medium ■ Cosmology: large-scale 
structure of Universe 



1 Preamble 

The present chapter provides the general framework of the physics of clusters of galaxies: it 
outlines how clusters connect to several astrophysical issues, from cosmology to the forma- 
tion of galaxies and stars. Some of these topics are discussed in the chapters of this volume; 
we refer to these chapters when appropriate. When the topic we mention here is not dealt 
with elsewhere in this volume, we refer to some of the most recent reviews. 



2 The observational framework 

2.1 Cluster components 
2.1.1 Mass content 

The discovery of clusters as concentrations of galaxies on the sky dates back to 1784, when 
Charles Messier mentions the Virgo clust er in his Connaissance des Temps (see the fasci- 
nating historical review by iBivianol 2000|). The first systematic optical surveys of clusters 
appeared almost two centuries later dAbelllll958l : IZwickv et al.lll961-68h . 

In 1930's, Zwicky, by assuming that the cluster in the Coma constellation, containing 
hundreds of bright galaxies (Fig.[T]iLjis in virial equilibrium, found that the mass required to 
bind the system gravitationally sh ould be roughly hundred times larger than the sum of the 



masses of the individual galaxies ( Zwicky| [l933l . 19371) . In fact, if the cluster is an isolated, 
spherically symmetric system in dynamical equilibrium, the virial theorem yields the total 
mass 

» 

where G is the gravitational constant, a is the dispersion of the galaxy velocities along the 
line of sight, and R is the cluster size 

where is the projected separation between galaxies i and j, and N is the number of 
galaxies; equation assumes that all the galaxies have the same mass. Clusters are not 
isolated syste ms and the correct application of the viri al theorem requires the inclusion of 
a surface term ( iThe & Whitdl 19861 : iGirardi et al.lll998h . However, correcting for this factor, 



See Biviano ( 1998) for a historical perspective on the Coma cluster. 
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Coma Cluster 

2.0 keV 




Fig. 1 Left panel: Optical image of the central part of the Coma cluster (23 x 23 arcmin) centered on the 
galaxy NGC4889; the bright galaxy on the right is NGC4874. Right panel: X-ray image (ROSAT) of Coma 
(164 x 150 arcmin) (from http://chandra.harvard.edu/photo/2002/0150/more.html I. 



and correcting also for departures from spherical symmetry, for galaxies of different masses 
and for the presence of interlopers, namely galaxies which are not cluster members but 
appear projected onto the cluster field of view, is not sufficient to fill in the discrepancy 
between the virial mass and the sum of the galaxy masses. We can reach the s ame conclusio n 
with a more sophisticated dynamical analysis based on the Jeans equations teiviandl2006h . 
Zwicky's result was the first indication of the existence of dark matter. 

With the advent of astrophysics from space in 1960's, the Uhuru s atellite revealed that 
clusters are the most luminous extended X-ray sources on the sky dGurskv et al.lll972l 
Fig. El. The origin of the X-ray luminosity L x was very early interpre t ed as thermal Brems - 
strahlung emis sion from a hot intr acluster plasma dFelten et alJI 19661 : iBvkov et alj|2008bl - 
Chapter 8, and lKaastra et al ] |2008l - Chapter 9, this volume): 



| W e(r)n io ns(r)A[r(r)]d 3 r : 



(3) 



where n e and «i ons are the electron and ion number densities in the intracluster medium 
(ICM) and A (T) is a cooling function. For temperatures kT >2 keV (k is the Boltzmann 
constant), when the ICM is almost fully ionised, we have A(T) °c T 1 / 2 . The existence of 
a hot diffuse X-ray emitting gas implies the presence of a deep gravitational potential well 
that keeps the gas confined. This evidence made unten able the hypothesis th at clusters are 
chance superpositions of galaxies or expanding systems i Ambartsumianl 1 96 lh . By assuming 
hydrostatic equilibrium and spherical symmetry, the cumulative mass within radius r is 



M(< r) 
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where \i is the mean molecular weight, m p the proton mass and p gas the gas mass density. 
This type of analysis confirmed that at least 70% of the cluster mass is dark llCowie et al.l 
Il987l) . and the rest is almost all in the diffuse ICM component: galaxies only contribute a 
few percent to the total cluster mass, as inferred by Zwicky in 1937. 




In this same 1937 paper, Zwicky also realized that one way to avoid the assumption of 
dynamical equilibrium when estimating the cluster mass is to use the deflection of light of 
background sources determined by the cluster gravitational potential well (Fig. for an 
axially symmetric cluster, the total mass projected along the line of sight within the distance 
r of closest approach of the light rays to the cluster centre is 

rc 2 

M(<r) = —a, (5) 

where c is the speed of light, and a is the deflection angle (e.g. ISchneiderj[2003) . The 
first indication of a gravit ational lensing effect in a galaxy cluster ap peared fifty years late r 
dLvnds & Petrosiar] 1986). To date, both strong gravitational lensing llHennawi et al] |2006). 



which creates multiple images of background sources, and weak gravitational lensing dDahld 
l2007h . which produces small induced ellipticity in the shape of background galaxies, indi- 
cate that, if General Relativity is correct on the scales of clusters, ~ 80 — 90% of the cluster 
mass is dark. 

An alternative method to estimate the cluster mass, both in their virial region and their 
outskirts, whic h does not rely on t he viri al eq uilibrium assum ption, is the caustic technique 
introduced by iDiaferio & Gelled d 19971) and iDiaferid (l999). This technique exploits the 
distribution of galaxies in redshift space to infer the galaxy escape velocity from the grav- 
itational potential well of the cluster, and hence its mass. In fact, in the redshift diagram, 
namely line-of-sight velocity vs. projected distance r from the cluster centre, cluster galaxies 
populate a trumpet-shaped region limited by two curves, named caustics, whose separation 
(r) at each radius r is proportional to the escape velocity. For a spherically symmetric 
cluster, the mass within radius r is 



M(<r) = ^V 2 (*)dx. 



(6) 
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Fig. 3 Contours of the SZ effect measurements with SZA overlaid on the X-ray images (XMM EPIC/MOS) 
of the clusters C1J1415.1+3612 (z= 1.03), C1J 1429. 0+424 1 (z = 0.92) and C1J122 6.9+3332 (z = 0.89). The 
filled ellipses show the FWHM of the SZ observations. See Muchovei et al. 1 2007) for details. 

All the mass estimation methods used to date, and described here, indicate that the 
dark matter contributes ~ 80 — 90% of the total cluster mass, the ICM contributes ~ 10 — 
20%, and the galaxies contribute less than a few percent. With the advance of X-ray spec- 
troscopy, it has beco me clear that the ICM is chemically enriched to ~ 0.5 Solar abundances 
jWerner et al.ll2008l - Chapter 16, this volume): given the relative mass contribution of the 
cluster components, it follows that the total mass in metals in the ICM is larger than the sum 
of the mass in metals of the individual galaxies. 

2.7.2 Clusters as probes for cosmology 

Right after the detection of the X-ray emission from clusters, ISunyaev & Zeldovichl dl972h 
realised that inverse Compton scattering of the Cosmic Microwave Background (CMB) pho- 
tons by the free electrons of the ICM can produce secondary anisotropies in the CMB maps: 
^/cmb/^cmb = f(v)y, where /(v) is a spectral function and 



is the Comptonization parameter; here, m e is the electron mass, <7t the Thomson cross- 
section, and the integral is along the line of sight I. Typical ICM temperatures kT and elec- 
tron number densities « e in clusters yield 4/cmb /^cmb ~ 10~ 4 (Fig. [3]>. Only over the last 
few years the sensitivity of radio telescopes and the control of systematics have been suf- 
ficient to measur e reliably this Sunyaev- Zeldovich (SZ) effect and its spectral signature in 
massive clusters dBonamente e t al. 2006). The major advantage of the SZ effect is that the 
amplitude of the CMB fluctuation is independent of the cluster redshift z, because the CMB 
intensity increases as (1 + z) 4 and compensates the (1 +z)~ 4 decrease of the SZ "surface 
brightness". This pr operty has recently bo osted an entire research field of cluster surveys 
with the SZ effect dCarlstrom et al.1 120020. In fact, "flux-limited" SZ cluster surveys are 
almost mass-limited cluster surveys dBarflettl 120061) . This property is extremely relevant, 
because the mass of clusters is a sensible cosmological probe, as we explain below. 

If clusters are in virial equilibrium, we can derive simple relations between their global 
properties, namely mass, galaxy velocity dispersion, number of galaxies (richness), X-ray 
luminosity, ICM temperature, and so on. By considering for example the virial relation 




(7) 
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3kT I (2[im p ) = GM/R, the scaling relation between the total mass M and the gas tempera- 
ture kT reads 

kT = 3.229 (±)(±) lft ( 14 M V /3 keV (8) 

where 8 is the average cluster overdensity w.r.t. the critical density 3Hq j (8kG) of the uni- 
verse, with Hq = lOOh km s _1 Mpc~' being the Hubble constant at the present time. Quan- 
tities in equation |8]l are normalized to typical cluster values. Analogously, we can write the 
total X-ray luminosity 

* (iPPis)'"' 1 ""' • (9 > 

where we approximated the cooling function A(T) = 0.843 x 10~ 23 (W/keV) 0,4 erg cm 3 
s _I at kT > 1 keV, which holds for gas with poor metallicity, and assumed n e = «[ ons = n = 
/ gas p/( J Um p ), where / gas is the fraction of the cluster total mass in the ICM and p is the 
cluster total mass density. 

These scaling relations provide a straightforward method to estimate the cluster mass 
from a directly observable quantity, like the X-ray temperature. Quantities related to X-ray 
observations have received particular attentio n because of their r obustness and the relatively 
simple physics involved in their correlations jRosati et alj|2002h . 
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Cluster masses are extremely relevant to constrain the cosmological model, because 
clusters populate the exponential tail of the mass function of galaxy systems. If the power 
spectrum of the primordial perturbations of the mass density field in the early universe is 
a power law with index n, the number per u nit volume \dn(M)/dM]d M of galaxy systems 
with total mass in the range (M, M + dM) is dPress & Schechteil[l974h 



dM JkM 2 V 3/ \mJ f 



(«+3)/6 



M \ (" +3 >/ 3 



dM, (10) 



where p is the (constant) comoving mean mass density of the universe and M* is a param- 
eter depending on the normalisation of the power spectrum o% and on the structure growth 
factor, which, in turn, depends on time, the cosmological density parameter Q m and the 
cosmological constant Q.a - M* increases with time when n > —3 and M» ~ 10 14 A — 1 Mq 
at the present epoch. Clearly, since the exponential cut off dominates the mass function at 
M > M*, the evolution of the cluster number density is a very sensitive indicator of the power 
spectrum normalisation and of the cosmological parameters (Fig. |4j). In practice, the appli- 
cation of this idea requires moder n versions of the Press- Schechter mass function which are 
more sophisticated than eq. d 1 b ( Sheth & Tormen|[l999l) : moreover, care must be taken to 
deal with the degeneracy among the cosmological parameters (e.g. IVoitll2005r) . especially 
between Q m and Os, which follow a rela tion of the form cig = AQ~ a , with A w 0.5 and 
a « 0.5 (e.g. lReiprich & B ohrineer 2002). Nevertheless, by simply assuming that the mean 
mass density of the universe must be smaller than the mean density within clusters and larger 
than the mean density computed including only the richest clusters, Abell already in 1965 
estimated £l m ~ 0-2 , a value remarkably close to the currently accepted estimate Q m = 0.26 
JSpergel et atll2007l) . 

In most mass estimation techniques, the assumption of dynamical equilibrium is fun- 
damental for obtaining accurate estimates. There is however overwhelming evidence that 
many clusters are out of equilibrium, as described in the following sections. 



2.2 Evidence of non-equilibrium and formation processes 
2.2.7 Cluster galaxies 

Elliptical and lenticulars (early-type) galaxies tend to reside in high-density regions, whereas 
spiral (late- type) galaxies are more common in low-density regions. This morphology-den- 
sity relation has been known since a publication of Shapley in 1926, and possibly earlier. 
The first statistical analysis based on a large sample of 55 clusters showed that, in the local 
universe, the fraction of spirals is 8 0%, 60% and 0% in the field, outskirts and cores of 
clusters, respectively jDresslej ri980). Subsequent work confirmed the connection between 
galaxy properties and local density: spirals in clusters are H I- and dust-deficient, have larger 
metallicity and show an excess of radio continuum emission (see iBoselli & Gavazzill2006l . 
for a review). 

At higher redshifts (z ~ 1), the fraction of spirals a nd irregulars in high-density re- 
gions increases at the expense of the fraction of lenticulars dPostman et al.l20 05; Smith etall 
|2005[) . Larger gas reservoirs and galaxy close encounters and collisions, which are more fre- 
quent in dense environments and favour the transformation of galaxy morphology, increase 
the star formation activity, as it occurs in starburst galaxies; indeed mid-infrared observa- 
tions of luminous infrared galaxies show a more intense star formation activity in cluster 



8 




0.0 B , , i , , Fl 

0.1 1.0 

Fig. 5 Differential velocity dispersion profiles of galaxies brighter than Mr = —20.5 in the CNOC1 clusters 
(0. 18 < z < 0.55). Top panel: Profiles of blue (open squares) and red (filled squares) galaxies; the rest-frame 
colour threshold is B — V = 0.85. The dashed and the solid lines show the corresponding profiles of the 
blue and red galaxies in mock clusters simulated with the N-body+semi-analytic technique (Diaferio et al. 
2001; see Sect.fJ). Bottom panel: Velocity dispersion profiles of disk-dominated galaxies (open squares) and 
bulge-dominated galaxies (filled squares). See Diaferio et al. (2001) for details. 

environment at increasing redshift (Met calfe et alJl2005l : lBai et al]|2007b . Although, on av- 
erage, the properties of galaxies vary smoothly from the cluste r centre to the outs kirts, before 
reaching the average properties of galaxies in the field (e.g. iRines et alj|2005b . the galaxy 
properties within individual clusters can show significant spatial variati ons, like the differe nt 
slopes of the galaxy luminosity function in different regions of Coma jAdami et al.ll2007h . 

In clusters with irregular morphology, the major axes of the bright galaxi es tend to be 
aligne d both with eac h other and with the major axis of the parent cluster jPlionis et al.l 
2003), as a century ago lWolj dl902h already pointed out for the galaxies in Coma. An align- 
ment between the major axis of the bright central galaxy (BC G) or the cD galaxy and th e ma- 
jor axis of the parent cluster is also present, as first noticed bv lCarter & Metcalfg dl980h . Any 
alignment present in the initial conditions is expected to be erased by multiple galaxy en- 
counters occurring in the cluster environment. Therefore, the presence of alignment should 
indicate a young dynamical state of the cluster. 

Blue cluster galaxies show larger velocity dispersions than red galaxies (Fig.[5l[Biviano et all 

1 1997b . This difference is attributed to more elongated orbits dBiviano & Katgertll2004b . in- 
dicating that, unlike the red galaxies which are likely to be in dynamical equilibrium within 
the cluster gravitational potential well, blue galaxies are, on average, possibly falling onto 
the cluster for the first time. 

Intracluster star s detected with planetary n ebulae or as diffus e light in clusters at low 
jGerhard et aiT 2005) and intermediate redshift dZibetti et al.ll2005b can be the relics of tidal 
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Fig. 6 Left panel: Temperatur e map of the cluste r RX J 1347.5-1 145. The superposed contours show the X-ray 
surface brightness. See Gitti & Schindler 1 2004) for details. Right panel: Map of the me tallicity abundance of 
the central region of A426, the Perseus cluster; brighter regions are more metal rich. See Sanders et al. ( 2005) 
for details. 



interactions suffered by falling galaxies dCovone et al. 2006) and/or the relics of merg- 
ing of galaxies alrea dy in place in the forming g ravitational potent ial well of the cluster 
(Mur anteet alll2007l) . as observed in real clusters ( IRines et al.ll20 07); intrac luster stars can 
contain from 10 to 50 % of the mass in stars in the core of clusters and can show distinct 
kinematic substructure dArnaboldi et"alll2006l) . 



2.2.2 Substructure, mass accretion and non-thermal phenomena 

A large fraction of clusters shows the presence of substructure both in their galaxy distri- 
bution and in their X-ray emission morphology. The fraction of clusters with substructure 
depends on the cluster sample and on the substructure identifi cation technique, but it is sub- 
stantial, lying in the range 30-80% (see iRamella et al. 2007, and references therein). The 
X-ray band also shows patchy temperature teelsole et al.ll2005j) and metallic ity maps (e.g. 
lHavakawa et"ai1l2006l : iFinoguenov e t al. 2006, Fig. HO see I Werner etaill2008l - Chapter 16, 
this v olume). Moreover, t he X-ray morphology is increasingly irregular with increasing red- 
shift jjeltema et ai]|2005l) . Where the angular resolution is low and detailed maps of the X- 
ray surface brightness cannot be obtained, a clear indication that, in some cl usters, the ICM 
migh t be out of equilibrium comes from the hard, possibly thermal, excess (Rephaeli et al. 
120081 - Chapter 5, this volume), which appears when the X-ray spectrum is fitted assuming 
a plasma with single temperature and metallicity. In these cases, two or more temperatures 
are required to yield a reasonable fit. 

All these pieces of evidence suggest that clusters accrete matt er from the surrounding 
regions. The cluste r 1ES 0657-55.8, dubbed "the bullet cluster" drvlarkevitch et al.l 12002 : 
Barr ena et alj 120021) . is a spectacular example of a merging cluster at intermediate red- 
shift (z = 0.296). At low redshifts, there are many other examples of merging clusters or 
galaxy groups falling onto clusters along filamentary structure. To mention a few: Coma 



10 



dColless & DunrJI 19961: 1 Adami et alJEoOsh . A 521 dFerrari et al.ll2003h . A 754 dHenrv et alJ 
l2004h. A 2 199 dRines et alJl200lh . A 22 19 dBoschin et alJl2004h . A 3560 in the Shapley su- 
percluster dBardelli et alj2002h . A 3921 dFerrari et al.l2005h. At high redshift the numbe r of 
irregular and merging clusters increases substantially dRosati 2004: lKodama et al]|2006h . 

Accretion of matter is a natural explanation for the patchy morphology of the X-ray 

emiss ion and metallicity m aps; it can also cause most of the turbulence in the ICM dSchuecker et al.l 
I2004T) and generate shocks (Bykovet al]|2008cj - Chapter 7, this volume), which can be effi- 
cient generators of high-energy cosmic rays. 

Non-thermal processes in the ICM appear in the X-ray and radio bands ( Repha eTi et alJ 
2008 - Chapter 5 , this volume). Non-thermal hard X-ray emission is indeed observed in a 
few systems (e.g. lPetrosian et alJl2006l ; lFusco-Femiano et alj|2007 | . and references therein) . 
However, the existence of both the hard X-ray emission excess liossetti & Molendill2004h 
and the soft X-ray and UV emission excesses, regardless of their thermal or non-thermal 
origin, are still deb ated (e.g. llieu & Mittazl2005l : lNevalainen et alj2007l : IWerner et al.l2007l : 
iDurret et alJl2008L - Chapter 4, this volume). 

At radio frequencies, the cluster emission can be both extended and associated with 
galaxies. One can infer the direction of the bulk motion of gas clouds within the ICM from 
the deviations that the I CM kinetic pres sure exerts on the radio jets emerging from active 
galactic nuclei (AGN) dBurns et ai]|2002h . 

The extended radio emission has a steep spectrum and can appear with two different 
morphologies: either the emission comes from the cluster centre (halo), is ex tended and reg- 
ular, or the emission com es from the cluster outskirts and is elongated (relic) ( Fere tti & Giovanninil 
l2007l : lFerrari et al.l2008l - Chapter 6, this volume). By combined radio and hard X-ray obser- 
vations or by Faraday rotation measurements one can infer the presence of magnetic fields in 
the ICM with inte nsity ranging from a few tenths to several jx Gauss (e.g. iBlasi et al]|2007l ; 
iFerrari et al . 2008, - Chapter 6, this volume). 

Both the radio and the hard X-ray emissions originate from relativistic electrons: the 
former is synchrotron r adiation, the latter is emission due to inverse Compton scattering 
with the CMB photons IPetrosian et al.1120081 - Chapter 10, this volume ]. Among the many 
accel eration mechanisms proposed to produce the relativistic electrons IPetrosian & Bvkovl 
2008 - Chapter 11, this volume], shocks due to the accretion of matter from the surrounding 
regions seem to be one of the most efficient. 



2.3 Connection with the large-scale structure 
2.3.1 A morphology -density relation for clusters 

It has been known for a long time that clusters are not randomly distri buted in the uni- 
verse but form large concentrations named superclusters (Shaple"vl ll933l) . The connection 
between clusters and the lar ge-scale structure has bec ome evident with extended ga laxy red- 



shift surveys, like the C fA dde Lapparent et al.ll 19861) . the SDSS dYork et alj|2000l) , and the 



2dF dColless et ai]|200lh redshift surveys. 

To quantify the distribution of objects in space, the simplest quantity to compute is the 
two-point correlation function |(r), which yields the probability of finding two objects at 
separation r in excess to a Poisson distrib ution. Clusters correlate according to the correla- 
tion function £, = (r/ro) 7 , with y w — 1.8 dBahcall & Soneiralll983h . The correlation length 
ro is proportional to the cluster richness: richer clusters are rarer and their mean interclus- 
ter separation d larger. One finds ro °c d 5 , and ro w 20 — 25/j -1 Mpc for rich clusters. 
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Fig. 7 Two-point correlation function of the 903 APM cluster sample (open circles) and of a subsample of 
clusters with significant substructur e (filled circ les). Insert: Correlation length rg vs. the significance level o~ 
of substructure. See Plionis & Basilakos (2002) for details. 



However, this ro — d relation cannot be fully responsible for the large range, from ~ 10 to 
~ 30/t Mpc, covered by t he correlation length rp wh en one considers clusters with in- 
creasing substructure (Fig. [7] Plionis & Basilakos 2002). Moreover, cluster shapes are gen- 
erally elongated, as first quantified bv lSastrvl (1968), and close clusters have their major axes 
aligned (Binggel fl982h . In the APM cluster sample, containing 903 objects, this alignment 
is stronger when the fraction of substructure is larger dPlionis & Basilakos 2002). 

This connection between the properties of clusters and their lar ge-scale environment is 
also apparent in the X-ray band: the X-ray luminosity o f clusters (Bohring er et al. IHqH) 
and the irregularity of the X-ray surface brightness maps dSchuecker et alfeOOll) are larger 
in environments with a higher cluster nu mber densit y. 

These features appear in a scenario ( Wes nl 19941) wher e clusters form by episodic mass 
accretion along preferential directions dColbergetal.ll999h . in agreement with the currently 
accepted model of the formation of cosmic structure (see Sect.[3]l. Future estimates of clus- 
ter peculiar velocitie s based on the kinetic SZ effect, when corrected of their systematics 
(Diaf erio et al.l f2005). will further test whether the correlation between the evolution of 
the peculiar velocity and the local density predicted by this model dPiaferio etail feoOO; 
ISheth & DiaferiofeoOlb is correct. 



2.3.2 The surrounding WHIM 

The optical dept h of the intergalac tic medium at z > 2 derived from the Lyman- a forest 
of QSO spectra jRauch et al] 1 1997b . the abundances of light elements combined with the 
predic tions of the standar d Big Bang nucleosynthesis model and the power spectrum of the 
CMB dSpergel et alj2007b are all consistent with a density of baryonic matter Q\,h 2 ~ 0.022. 
At low redshift, baryonic matter within galaxies and with in the diffuse hot med ium of galaxy 
clusters and groups only accounts for ~ 15% of this Q.\, (F ukugita et al.|[l993) . 

The hierarchical clustering model solves this discrepancy by predicting that 50% of the 
baryons permeating the matter distribution is in a warm phase with kT in the range 0.01 — 1 
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keV, overdensity in the range 0. 1 — 10 and median ~ 10 — 20, and m edian metallici ty ~ .2 
Solar, with a large spread, for oxygen dCen & Ostrikeril 19991 . 120061 : | Berto ne et al ] |2008l . - 
Chapter 14, this volume). The WHIM overdensity increases with increasing dark matter 
overdensity. Therefore, the outskirts of clusters and their environments are expected to be 
sources of a thermal soft X-ray emission. However, both the expected low surface bright- 
ness of the WHIM and the intergalactic and Galactic absorption make this observation very 
challenging. Neve rtheless, there are claims of detection in many clusters, incl uding Coma 
(iTakei et al.ll2007l) . In fact, rather than detecting the thermal X-ray emission l lDurret et al.l 
120081 - Chapter 4, this volume) directly, the most promising technique is searching for ther- 
mal UV and X-ray absorption lines in the spectra of bright background quasars or blazars 
(Richter et al. 2008 - Chapter 3, this volume). Even with this strategy, current instruments 
(Paer els et al.l20 08 - Chapter 1 9, this volume) must be used at their sensitivity limits and the 
most plausible detections (e.g. iNicastro et al.ll2005h are indeed still debated ( iKaastra et alj 
l2006h . Firm detections are however of great relevance, because the existence of the WHIM 
and the estimate of its properties would be an extraordinary confirmation of the current 
model of structure formation. 



3 The modeling framework 

3.1 Players: gravity and beyond 

Modeling the formation and evolution of galaxy clusters requires the set-up of a complete 
cosmological c ontext. According to the current paradigm of structure formation (as re- 
viewed, e.g., in ISpringel et alJl200a) . inflation amplifies quantum fluctuations, which are 
present in the dark matter density field of the early universe, to cosmic scales; at the time 
of matter and radiation density equality, these fluctuations start growing by gravitational 
instability. After radiation and ordinary baryonic matter decouple, radiation pressure stops 
supporting the baryonic density perturbations against self-gravity and against the pull of the 
dark matter gravitational potential wells that have already formed: the dark matter inhomo- 
geneities thus accelerate the collapse of ordinary matter. 

The pre-inflationary Gaussian quantum fluctuations generate a scale-free power spec- 
trum of the primordial density perturbations. At later times, the dark matter particles, which 
have sufficient momenta to stream out of the denser regions, set in a characteristic scale in 
the power spectrum, which is proportional to the typical velocity of the dark matter par- 
ticles themselves; below this scale the density perturba tions are damped. The existence of 



quasars ana galaxies at high reasmit (ranetal. 2UU4), tne power spectrum ot the LMB 
jSmootetal.ll 19921 ; Ide Bernardis et al. 2000; Hanany et al. 2000), and the power spectrum 
of the large-scale distribution of galaxies dTegmark et al.l2004l) all agree with the hypothesis 
that the dark matter is "cold", namely that the dark matter particles become non-relativistic 
at very early times and the damping scale is <C Ih^ 1 Mpc. In these Cold Dark Matter (CDM) 
models, larger systems form by the aggregation of smaller systems and galaxy clusters are 
the nodes of the network of sheets and filaments which constitutes the matter density field 
(Fig. [8}. Clusters can form at reasonably recent times, with the actual formation time de- 
pending on the cosmological parameters. 

In CDM models, the scale-free properties of the primordial power spectrum are pre- 
served on scales larger than the scales of galaxies. This feature implies that the dark matter 
halos of galaxies, groups and clusters evolve self-similarly. If gravitation is the only force 
playing, the physics of the ICM is ruled by adiabatic compression and shock heating, and we 
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Fig. 8 Map density at z = of the gas in a slice of the N-body/hydrodynamics simulation of Borgani et al. 
(2004). The volume shown is 192/;~' Mpc on a side and 1 2/i -1 Mpc thick. The high-density regions are 
white, the low-density regions dark. See lBorgani et al l 1 2004) for details. 



can derive simple relations between cluster properties at different redshifts jKaiser| [T986). 
The cluster virial mass scales as M v ;,- °c p c (z)A c (z)R 3 ', where R is the cluster size in physical 
units, p c (z) = 3// 2 (z)/(8?rG) is the critical density of the universe, and A c (z) is the cluster 
density in units of p c (z);Elpc(z) scales with redshift z as p c (z) « E 2 {z) = X2 m ( 1 + z) 3 + (1 — 
Q m - Q A )(l +z) 2 + C2 A - The cluster size thus scales as R °c mV 3 4 l/3 (z)£" 2 / 3 (z), and 
the temperature asr« M vil /R °c mV 3 4 c 1/3 (z)£ ,2 / 3 (z); thus, for example, one finds for the 
X-ray luminosity due to Bremsstrahlung emission: 

' n e n lons A(T)dV oc M 2 m R- 3 T^ 2 oc ^ E 1 ^ (z)^ (z) (12) 

or for the Comptonization parameter 

y oc [ n e Tdl oc M vir R- 2 T oc M vir E 2 {z)A c (z) . (13) 



These scaling relations only partially agree with observations; in fact, (1) the relative 
contribution of non-gravitational sources (AGN, galactic winds, supernovae, turbulence) to 
the ICM energy budget depends on the mass of the cluster and affects the slope of the scaling 
relations; and (2) the merging history of the cluster introduces an intrinsic scatter in the 
relations dCavaliere et al.|[l998l) . More over, clusters at h igh redshift have scaling relations 
steeper than the self-similar prediction tettori et alj|2004l) . 



2 A well-known approximation is 



where w = Q m (z) - 1 iBrvan & Normanll998l) . 
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The deviations from the self-similar scaling relations are clearly related to non-gravi- 
tational processes. The interplay between gravitational and non-gravitational processes is 
highly non-linear and treating it appropriately requires sophisticated modeling. Over the last 
thirty years two main approaches have been followed: numerical simulations dDolag et alj 
l2008al - Chapter 12, this volume) and an analytic (or semi-analytic) approach. Over the last 
ten years, the two approaches have been combined. 



3.2 Making clusters: recipes, ingredients and tools 

For the last thirty years, numerical simulations have represented the main tool of theoretical 
astrophysicists to investigate the non-linear evolution of cosmic structure and compare the 
theoretical results with observations. The first comput ational device fo r simulating galaxy 
interactions used an ingenious system of light-bulbs jHolmberdfl94lh . Numerical simu- 
lations of astrophysical systems have now reached tremendous sophistication dPolag et alj 
2008a - Chapter 12, this volume) when compared to the first numerical integration of the 
equations of motion of a few stars in a star cluster dvon Hoerne3 [T960). Current numeri- 
cal models reproduce the gros s features of the larg e-scale struc ture, and the distrib ution of 
galaxies and diffuse baryons (Borg anietal.ll2008al - Chapter 13: lBertone et alj2008l - Chap- 



ter 14 , this volume), including the statistical properties of the Lya forest I Springel et alj 
l2006h . 

Until the 1980's, numerical simulations relied on iV-body techniques only to follow 
the evolution of the dissipationless component of cosmic structure. Particles in the sim- 
ulations were first identifi ed wi th galaxies and later on with dark matter particle tracers. 
Peeble] dl970h and lWhitd dl976f) were the first to simulate the formation of a galaxy cluster 
with vacuum boundary conditions. To obtain a realistic distribution of matter, however, it 
is necessary to simulate a cosmological volume, with periodic boundary conditions, which 
is large enough to be representative of the universe. The first simulation of this kind, with 
a cold dark matter power spectrum of the initial density perturbations, demonstrated that 
dark matter must be c onstituted of cold collisionless particles, rather than massive neutri- 
nos dPavis et"al]|l985 ). Recent dark matter only simu lations yield clusters whose number 



density evolution fevrard et alj|2002h . virial relatio ns fevrard et al.1 2007), and connection 
between their shape and the large-scale structure dAltav et alj|2006l : iLee & Evrardll2007h 

resemble the observations and can be used as probes of the co smological model. 

In the late 1980s, smoot hed particle hyd rodynamics (SPH) fevrard! 198 8: Hernq uist & Katzl 

1 19891) and mesh techniques dCen et al.l l 1990) were introduced to follow the evolution of the 
dissipative component. The presence of dissipative gas substantially increases the number 
of physical processes to be taken into account and therefore increases the complexity of the 
simulations. In fact, a priori, it is not guaranteed that the different techniques yield similar 
results and some code comparisons have been performed over the years to check the consis- 
tency of the various integration schemes ( IKang et alll994l : lFrenk et alJll999l : l0'Shea et alj 
120051) . 

The simpl est approxi mation is to consider the gas non-radiative and perform adiabatic 
simulations fevrardlll99Cl) . This model is able to reproduce the general trend of the self- 
similar scaling laws, but disagrees with observations in detail, yiel ding, in general, shal - 
lower X-ray luminosity-temperature and mass-temperature relations dNavarro et ai][T9 95b). 
To reproduce the observa tions more closely, and expla i n the deviations from the self-similar 
model, radiative cooling dThomas & Couchmanll 19921 : iKatz & Whitel fl993) and non-gravi- 
tational heating processes dNavarro et al.ll 1995bh must be included. More realistic models, 
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attempting to reproduce the thermal properties of the ICM and the stellar content of clus- 
ters, have to include explicitly all those intimately coupled processes which rule galaxy 
formation and th e thermodynamic s of the ICM: star fo rmation, energy feedback from super- 
novae explosion Jvaldarnini 2 0021) and ga lactic winds dMetzler & Evrardlll994l ; lDolag et alj 
2008a - Chapter 12: lBorgani et al. || 2008al - Chap t er 13, this volume), besides metal produc- 
tion and chemical en richment I IValdarninilEooilJSchindler & Diaferiol [20081 - Chapter 17; 
iBorgani et all2008bl - Chapter 18, this volume), which strongly affects the gas cooling rate; 
moreover, one needs to include magn etic fields and its associated non-thermal processes 
(Dolag et al.lll999l ; lDolag et alj J2008bh - Chapter 15, this volume). 

The galaxy population in the simulated clusters have shown some rough similarities with 
observations since the very first numerical experiments dEvrard et al.ll 19941) ; many other ob- 
served properties are reproduced at d ifferent levels of agreement, including the soft X-ray 
thermal emission dCheng et al ] l2005h . Nevertheless, in no simulation today the thermody- 
namical properties of the ICM, the stellar mass fraction and the galaxy luminosity function 
simultaneously match the observations. One of the most serious discrepancies with obser- 
vations is the overproduction of stars: the fr action of baryons in stars usually is ~ 20%, 
at least a factor of two larger than observed jBorgani et ai]|2004l) . A possible solution has 
been identified in the inhibition of the ICM cooling by the energy injection of AGN jets, as 
suggested by the X-ray cavities in correspondence of the radio emission ob served around 
cluster central galaxies (Mc Namara et al.ll2006l : iMcNamara & Nulsenll2007h : the presence 
of these cavities correlates with the X-ray temperature d rops in the centre of some clusters, 
early interpreted as cooling flows dDunn & Fab ian 2006); however, the details of the heating 
mechanism by AGN jets remain uncertain I VoifeOOSo T 

The main difficulty when treating most of the dissipative processes resides in the fact 
that they occur on scales much smaller than the spatial resolution of the simulation and 
must be approximated by phenomenological recipes: gas cooling, star formation and stellar 
evolution, which occur at scales <C 1 pc, have effects on galactic and extragalactic scales 2> 1 
Mpc; thus the simulations should model a spatial range of at least six orders of magnitude. 
These large dynamic ranges can be covered by neither a single numerical experiment nor 
experiments which rely on the resimulation technique, where an individual cluster extracted 
from a large-scale cosmological simulation is resimula ted at a much higher resolution with 
appropriate boundary conditions dNavarro et allll995ah . 

The large number of physical processes and the large computational resources required 
to run an /V-body/hydrodynamics simulation for a realistic modeling of cluster and galaxy 
formation have suggested a complementary strategy to explore more efficiently the param- 
eter space: the combination of semi-analytic model s with merger trees of dark matter halos 
extrac ted from dissipationless N-body simulations dRoukema et al.lfT997l ; iKauffmann et alj 
1999). This approach has reproduced the evolution and p resent-day spatial an d kinematic 
distribution of galaxies in the universe and in clusters dSpringel et al. I l2005h . The major 
shortcoming of this strategy is that severe uncertainties in the predictions can be introduced 
by the crude simplification of the physical processes and especially of the thermodynamics 
of the ICM, usually assumed to be isothermal and in hydrostatic equilibrium within each 
dark matter halo. Some improvements can be introduce d by describing the IC M thermal 
properties in terms of entropy generation and distribution dMcCarthv et al.ll2007h , but much 
work remains to be done. 
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